Diffusion-reaction mechanisms of nitriding species in Si02 
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We study using first-principles total-energy calculations, diffusion-reaction processes involved in 
the thermal nitridation of Si02. We consider NO, NH, N2 and atomic N in different charge states 
as the nitriding species in Q-quartz. Our results show that none of neutral species react with the 
Si02 network remaining at interstitial sites. Therefore, they are likely to diffuse through the oxide, 
incorporating nitrogen at near-interface (Si/Si02) regions. Whereas, charged species are trapped 
by the network, nitriding bulk Si02. For the diffusing species, we find that NH and atomic N show 
increasing diffusivities with temperatures, whereas for NO and N2 they are relatively constant. This 
result agree well with the finding of higher N concentration at the Si/Si02 interface obtained by 
thermal processing of Si02 in NH3 as compared with those obtained in N2O. Finally, we discuss 
spin-dependent incorporation reaction mechanisms of NH and atomic N with the Si02 network. 

PACS numbers: 66.30.Ny, 71.55.Ht, 71.15.Nc, 71.15.Mb 



I. INTRODUCTION 

Nitrided silicon oxide or oxynitride is currently the 
near-term solution to substitute Si02 as the gate insula- 
tor material for the ultrathin metal-oxide-semiconductors 
(MOS) technology. Recent studies have suggested that 
the performance of oxynitridc-based MOS devices de- 
pends both on the depth concentration and on the dis- 
tribution of nitrogen into the gate oxide. ^ According to 
these studies, the best nitrogen profile for an ultrathin 
gate oxide would have: (i) a small nitrogen concentra- 
tion near the Si02/Si interface in order to reduce degra- 
dations by hot electrons,^ and (ii) a larger N concen- 
tration near the interface between the dielectric and the 
polycrystalline silicon (poly-Si) gate electrode, in order 
to minimize dopant diffusioui^ 

The growth of ultrathin oxynitride films depends 
strongly on the reactant agent {e.g., N2O, NO, NH3) 
and the technique used. Nitrogen can be incorporated 
into Si02 using either thermal oxidation and annealing 
or chemical and physical deposition methods. Thermal 
nitridation of Si02 in N2O generally results in a rela- 
tively low N concentration at near-interface (Si/Si02) 
region. The nitriding species are originated in the de- 
composition of the N2O molecule occurring at typical 
oxidation temperature, these species being the NO and 
N2 molecules4i^ On the other hand, nitrogen incorpo- 
ration in Si02 can also be performed via annealing in a 
NH3 atmosphere, resulting in a relatively high N concen- 
trations into the films. ^ This method provides both near- 
interface and near-surface nitridation, which suggests dif- 
ferent mechanisms for N incorporation or different nitrid- 
ing species derived from NH3. The above thermal pro- 
cessing are performed at high temperatures (> 800 °C). 
Plasma nitridation is a promising method for making ul- 
trathin oxynitride at lower temperatures (~ 300-400 °C). 
Higher N concentrations and controlled distribution can 
be attained with this method,^ typically using ions and 
radicals derived from N2 and NH3 as nitrogen sources 
Additionally, nitridation by energetic particles (N ions) 



provides high N concentrations much closer to the near- 
surface region, with little or no nitrogen at the Si02/Si 
interface.! Although the nitridation mechanisms to make 
ultrathin Si02 films are well known at the few layer level, 
less is known about the diffusing mechanism of the ni- 
triding species and their reactions with the oxide at the 
atomic level. 

In this work the energetics and diffusing properties 
of the Si02 thermal nitridation are studied from first- 
principles total-energy calculations. We have considered 
the N2, NH and NO molecules as well as atomic N in 
different charge states as the nitriding species reacting 
and diffusing through the Si02 network. The outline of 
the paper is as follows. Section II describes the theo- 
retical procedure. Our results for the interaction of the 
nitriding species with the Si02 network is presented in 
Sec. IIIA. The diffusivity of the nitriding species through 
Si02 is presented in Sec. IIIB. In Sec. IIIC we discuss the 
spin-dependent diffusion reactions of NH and atomic N 
with the Si02 network. Finally, in Sec. IV we show our 
conclusions. 



II. THEORETICAL METHOD 

Our calculations were performed in the framework of 
the density functional theory,^ using a basis set of nu- 
merical atomic orbitals as implemented in the SIESTA 
codeiiS We have used a split-valence double-C basis 
set plus the polarization functions as well as standard 
norm-conserving pseudopotentialsfii For the exchange- 
correlation potential we adopt the generalized gradient 
approximationii^ We used a 72-atom a-quartz supercell 
and the F point for the Brillouin zone sampling. The po- 
sitions of all the atoms in the supercell were relaxed until 
all the force components were smaller than 0.05 eV/A. 
We also consider neutral and singly charged species, 
where the neutrality of the cell is always maintained by 
introducing a compensating background charge. Spin- 
polarization effects are included throughout the calcula- 
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FIG. 1: Local equilibrium geometries for singly charged ni- 
triding species after reacting with the Si02 network, (a) neg- 
ative NO. (b) negative and (c) positive NH, respectively, (d) 
negative and (e) positive atomic N, respectively, (f) the per- 
fect Si02 network. 

tion since they are important for the correct description 
of atomic and molecular reaction processes in Si02 We 
initially study the energetics and structural properties of 
N2, NH, NO, and atomic N in the largest interstitial site 
of a-quartz for three different charge states 0, — ). We 
explore possible reactions that these species may undergo 
with the network as well as their diffusivities in a-quartz. 

III. RESULTS AND DISCUSSION 

A. Incorporation of nitriding species into the Si02 
network 

Thermal processing of Si02 in N2O and NO atmo- 
spheres shows that the NO molecule is the species which 
diffuses through the oxide, being incorporated at the 
Si/Si02 interface.® In our calculations we find that neu- 
tral NO and N0+ do not react with the Si02 network, re- 
maining at the interstitial position and being, in this way, 
able to diffuse through the oxide as experimentally ob- 
served. However, NO^ is trapped by the network forming 
a structure where the N atom of the molecule is bound 
to two fourfold coordinate Si atoms of Si02, as shown 
in Fig. ^a). The binding energy of N0~, calculated as 
the difference in energy between the bound and the free- 
interstitial configurations for the same charge state, is 



found to be 1.3 eV. The Si-N and N-0 bond lengths are 
2.26 and 1.23 A, respectively, whereas the Si-Si distance 
decreases by about 6% with respect to the equilibrium 
Si02 structure. For interstitial NO and N0+ the N-0 
bond lengths are 1.17 and 1.13 A, respectively. 

The NH radical seems to be the most likely diffusing 
species after thermal dissociation of NH3 according to 
the dissociative reactions: NH3 NH2 -I- H — 4.5 eV and 
NH2 ^ NH + H - 3.9 eV. Although NH2 may be another 
diffusing species for the interface nitridation, it is a rela- 
tively larger molecule and it may suffer additional disso- 
ciations at near-surface Si02 during thermal processing. 
According to our results, neutral NH in the gas phase 
has a spin triplet (S" =1) ground state whereas the singlet 
(5=0) states is 2.06 eV higher in energy. However, in- 
side the largest interstitial site of a-quartz, this difference 
in energy between both spin configurations decreases to 
1.04 eV, due to the interaction with the crystal field. We 
find that triplet NH does not react with the Si02 net- 
work remaining relatively inert at the interstitial sites. 
However, singlet NH reacts quickly with the oxide form- 
ing a structure where the N atom of NH binds to both Si 
and O atoms of the oxide. This incorporation reaction 
is highly exothermic with an energy gain of 2.04 eV. In 
Section HI C we discuss the possibility that NH might be 
incorporated in the Si02 network via a triplet-to-singlet 
spin exchange. 

Singly charged NH molecules are highly reactive in the 
Si02 network forming bound configurations as shown in 
Fig. ^b) and Fig. ^c). The binding energies for NH~ 
and NH+ are 3.2 and 2.7 eV, respectively, which sug- 
gest that these configurations are very stable favoring a 
Si02 near-surface nitridation. We observe that the NH^ 
bound configuration is similar to that previously found 
for N0~. However, the relaxation of the Si02 network is 
larger for the NH" bound configuration [see Fig. ^b)]. 
Here, both Si-N bond lengths are about 2.1 A, whereas 
the Si-Si distance decreases by about 14% with respect 
to the perfect Si02 showing a strong lattice relaxation. 
The N-H bond length is found to be 1.04 A. In the bound 
NH"*" configuration [Fig. QJc)] , the N atom binds to an O 
atom of the oxide keeping its bond to the H atom. The 
0-N and N-H bond lengths are 1.43 and 1.05 A, respec- 
tively, forming a 0-N-H angle of 103°. 

The N2 molecule is a product of the N2O gas decom- 
position at typical oxidation temperatures, being intro- 
duced in this way into Si02 in thermal processing.^ N2 
(and atomic N) may also be introduced at lower temper- 
atures by plasma assisted methods.^ Our results for N2 
in Si02 show that this molecule does not react with the 
oxide being relatively inert for the three charge states 
considered. Therefore, it may diffuse easily through the 
oxide reacting with the silicon at the Si/Si02 interface 
or escaping from the films. 

The N atom in free space has a quartet (iS'=3/2) 
ground state. The difference in energy with respect to 
the doublet (5' =1/2) state is calculated to be '^3 eV. 
However, for a N atom in the largest interstitial site of 
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a-quartz, this energy difference decreases to 0.76 eV. We 
find tliat quartet N does not react witfi the Si02 network 
suggesting that it would be a diffusing species. However, 
doublet N reacts with the oxide being incorporated into 
Si02 network forming a Si-N-0 bond. This reaction is 
exothermic with an energy gain of 0.76 eV. In Section III 
C, we discuss the possibility that atomic N might be in- 
corporated in the Si02 network via a quartet-to-doublet 
spin exchange. 

Singly charged N atoms are highly reactive in Si02. 
binds to a fourfold coordinate Si atom forming an 
additional Si-N bond of 1.85 A, as shown in Fig. ^d). 
The N atom has a binding energy of 2.2 eV. On the other 
hand, N"*" binds to an oxygen forming a 0-N bond, as 
shown in Fig.Ue). We find a 0-N bond length of 1.39 A 
with a N binding energy of 1.1 eV. 

As a general trend, we can say that because of the 
ionic character of Si02, charged species are trapped by 
the network. The positively (negatively) charged NO, 
NH and N species attach to O (Si) atoms of the oxide, 
forming strong bonds, with binding energies ranging from 
1 to 3 eV. This suggests a mechanism for the high-density 
bulk and near-surface nitridation. On the other hand, the 
neutral species in their ground-state spin configurations 
do not react with the Si02 network. Therefore, they are 
able to diffuse through the oxide. 
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FIG. 2: Total energy variations for the neutral nitriding 
species diffusing through a channel normal to the c-axis of 
Qf-quartz. (a) for the NO and N2 molecules, (b) for the NH 
molecule and atomic N. The arrows indicate the point where 
the nitriding species pass closer to the Si02 network and the 
zero distance is the center of the largest void in a-quartz. 



B. Diffusion of nitriding species in Si02 

We calculate the diffusion coefficient for the neutral 
species hopping through the larger channel of a-quartz, 
normal to the c-axis. To calculate the migration barriers, 
we have followed the same procedure as described in our 
previous calculation of the O2 diffusion in a-quartz We 
fix one atom of the molecules at several points along the 
pathway joining neighboring interstitial sites of a-quartz 
where the distance between them is ~ 5 A. All other 
atoms of the systems are allowed to relax. Our results 
for the total energy variations along this diffusion path 
are shown in Fig. |21 Here, the zero distance indicates 
approximately the center of the largest interstitial site 
and the arrow, the place of closest approach to the net- 
work. As a general trend, we find an anisotropic energy 
profile through the migration path which can be associ- 
ated to the asymmetry of neighboring interstitial sites. 
Fig. [5fa) shows the calculated total-energy variations of 
diffusing NO and N2 ; we observe that these species show 
similar global minima and height barriers along the path- 
way. However, N2 exhibits two local minima not found 
in the NO pathway, showing that the interactions be- 
tween the crystal field and the molecular orbitals of the 
diffusing species may affect the diffusion energy profile. 
Nevertheless, we must consider that uncertainties are in- 
troduced by fixing one atom of the molecules in order to 
map the diffusion energy profiles. For NH and atomic 
N we find lower-energy barriers and shorter distances be- 
tween minima, as compared with NO and N2 (see Fig.|2l, 



suggesting higher diffusivities for these species. 

We estimate the thermal diffusion coefficient or diffu- 
sivity for the nitriding species using the Arrhenius law, 
given by 

DiT)^DoeM-Ea/kT). (1) 

Here, the prefactor is defined by Do — Piy/G for a three- 
dimensional migration path, where I is the hopping dis- 
tance between minima and v is the attempt frequency, 
kT the Boltzmann constant times the temperature, and 
Ea is the highest energy barrier of the migration path. 
J/ is calculated from the energy profile of each species at 
the interstitial site (approximately the zero distance in 
Fig-IUj using the harmonic approximation, whereas Ea 
and I are obtained directly from the same figure. Tab. d 
lists the values of these quantities as well as our results 



TABLE I: Energy barriers (Ea), hopping distance between 
minima (l) and attempt frequencies (u) as obtained from the 
migration energy profiles of Fig. |21 -Do is the calculated pre- 
factor coefficient of the diffusivities for each nitriding species. 



Species 


£a(eV) 


l{A) 


!/(xl0^^sec~^) 


Do(xlO"'*cmVsec) 


NO 


0.95 


4.7 


2.26 


8.32 


N2 


0.95 


4.8 


2.38 


9.14 


NH 


0.48 


2.8 


3.67 


4.79 


N 


0.60 


3.6 


3.41 


7.37 
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FIG. 3: Calculated difTusivities (D) of the nitriding species in 
a-quartz as a function of typical plasma assisted and thermal 
processing temperatures. 



for the prefactor of each nitriding species. It is worth 
pointing out that we have considered a very specific mi- 
gration pathway through a-quartz. Other pathways may 
result in dilFerent migration barriers and diffusivities. In 
fact, recent calculations have shown major differences in 
the diffusivity of atomic hydrogen in different crystalline 
structures of Si02«i^^ However, we find that our ap- 
proach is valid to obtain informations on the relative dif- 
fusivities of the species under consideration. 

Figure 121 shows the diffusivities of the nitriding species 
in a-quartz as calculated from Eq. ^ for a range of tem- 
peratures typically used in thermal and plasma assisted 
methods. Here we observe that NH and atomic N have 
the higher diffusivities in a-quartz with increasing tem- 
peratures as compared with NO and N2, which suggest 
that they would be the best thermally-activated diffusing 
species in Si02 and the most efficient for Si/Si02 inter- 
face nitridation, as compared with NO and N2. Accord- 
ing to our results, the ratio between NH and NO diffusivi- 
ties at 800°C is estimated to be -Dnh/^no — 90. Assum- 
ing that NH and NO are the main diffusing species, our 
results agree well with the experimental finding of higher 
N concentrations at the interface obtained by thermal 
nitridation of Si02 in NH3 as compared with thermal 
nitridation in NOi^ 



C. Spin-dependent diffusion reactions of NH and N 
in Si02 

As mentioned above, NH in the triplet state does not 
react with the Si02 network remaining inert at intersti- 
tial sites. However, singlet NH reacts easily with the ox- 
ide being incorporated into the network. Because of this 
spin-dependent reaction of NH with the Si02 network, 
we study the possibility that a diffusing NH may suf- 
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FIG. 4: Triplet {S = 1) and singlet {S = 0) total-energy 
curves along the reaction pathway {d in the inset figure) for 
the NH molecule approaching the Si02 network. The zero 
distance is the equilibrium position of NH in the largest in- 
terstitial site of a-quartz where its corresponding total energy 
is chosen as the zero energy of the system. 

fer a triplet-to-singlet spin conversion while approaching 
the oxide, resulting in its incorporation into the network. 
For this, we compute the potential-energy surface (PES) 
along a pathway joining the equilibrium position of in- 
terstitial NH and a Si-0 bond of the Si02 network, for 
both triplet and singlet spin states, as shown in Fig0| 
The triplet PES in the figure is depicted along a di- 
rection perpendicular to the diffusion pathway shown in 
Fig. 121 however, both pathways intercept at the inter- 
stitial equilibrium configuration of NH (zero distance in 
Fig. 21) . The probability that the system changes from 
the triplet to the singlet PES may be estimated by the 
Landau-Zener theoryiii This theory has been recently 
used to describe spin effects in the adsorption of O2 in 
the Si(OOl) surfacei^ and the O2 incorporation reaction 
in the Si/Si02 interface. -'^^ Considering that the NH is ini- 
tially in the triplet diffusion PES, and is evolving towards 
the crossing region with a velocity v, the probability for a 
conversion to the singlet PES (Pts) may be approximated 
by 

Pt, = 2[l-exp(-l/Vto|Fi-F2|)], (2) 

where V is the triplet-to-singlet spin-orbit matrix ele- 
ment of NH {X^Y.-^ b^Y+) of 65.1 cm^\i9 Fi and 
F2 are the slopes of the two PES at the crossing point 
and h is the Planck's constant, v is estimated from the 
NH center-of-mass thermal velocity at 800 °C, a typi- 
cal annealing temperature. -Fi and F2 are obtained from 
the triplet and singlet curves at the crossing region (see 
Fig-O- Thus, the probability for a triplet-to-singlet con- 
version is found to be Pts = 9.3x10"''. This means that 
triplet NH in a single passage through the crossing point 
has a small probability to change to the singlet PES. 
However, as the system is thermally activated, many pas- 
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sages of NH through the crossing point will be performed. 
In order to quantify the spin-conversion events for NH 
molecules approaching the Si02 network at a given tem- 
perature, we estimate the ratio between the rate of NH 
following the triplet diffusion PES (Ft) and the rate of 
NH changing from triplet to singlet PES at the crossing 
point (Pts)- Assuming that the thermally-activated sys- 
tems follow the Arrhenius law, the reaction rates can be 
written as 



1.5 



and 



Tt - To eM-^Et/kT) 



(3) 



(4) 



Here, Pq is the frequency prefactor; AEt is the en- 
ergy barrier for the triplet diffusion PES obtained from 
Fig. Elb) {AEt =0.5 eV) and AEp is the energy differ- 
ence between the lowest-energy equilibrium configuration 
of the triplet NH and the crossing point, obtained from 
Fig. H (A-Bp = 0.1 eV). For a temperature of 800 °C, we 
obtain for our calculated Pts that Tt/Tts — 14, which 
indicates that NH will proceed preferentially along the 
triplet diffusion PES. 

As the incorporation of NH into the oxide due to the 
spin exchange is also likely to occur (an average of one in 
14 events), we study three possible thermally-activated 
diffusion-reaction mechanisms for an incorporated sin- 
glet NH at 800 °C: (i) NH diffuses along the Si02 net- 
work by a hopping mechanism. The rate is estimated 
by Eq. |3| substituting AEt by a calculated hopping en- 
ergy barrier {AEh = 2.2 eV). We find a hopping rate of 
F/i = Fo(4.6x 10~^^). (ii) NH is released from the network 
following the singlet PES. This reaction rate is estimated 
by 



Ts = ro[cxp{-AEs/kT)]x{l-Pts), 



(5) 



where AE^ is the singlet PES energy barrier ob- 
tained from Fig. H {AEs=2.0 eV). We find that 
Fs =ro(4.0xl0"i"). (iii) NH is released from the net- 
work changing from the singlet to the triplet PES at the 
crossing point. The reaction rate is estimated by Eq. ^ 
where AEp is now the difference between the energies 
at the singlet-triplet crossing point and at the incorpo- 
rated NH equilibrium configuration, obtained from Fig. 01 
(A£;p = 1.05 eV). We find that F^t = ro(l.l x 10"^). Ac- 
cording to these results, the most likely diffusion mech- 
anism for an incorporated singlet NH in Si02 would be 
mediated by the reaction (iii), i.e., NH leaves the net- 
work changing from the singlet to the triplet PES at the 
crossing point and then diffuses through Si02 as a triplet 
NH. The two other mechanisms are less likely, as can be 
verified by the ratios Fst/F/i ~ 240 and Fst/Fs ^ 28. 

Similarly to NH, we find a spin-dependent diffusion re- 
action for atomic N with the Si02 network. According 
to our results, the N atom in the ground quartet state 
does not react with the network remaining at intersti- 
tial sites, whereas, in the doublet state, the N atom is 
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FIG. 5: Doublet {S = 1/2) and quartet {S = 3/2) total- 
energy curves along the reaction pathway {d in the inset fig- 
ure) for the N atom approaching the Si02 network. The zero 
distance is the equihbrium position of the N atom in the 
largest interstitial site of a-quartz where its corresponding 
total energy is chosen as the zero energy of the system. 



incorporated between a Si-0 bond forming a structure 
similar to the peroxyl bridge of oxygens2&. Therefore, 
we examine the possibility that the N atom, initially in 
the quartet state, changes the spin configuration to the 
doublet one, resulting in its incorporation into the oxide. 
Following the same procedure described above for NH, 
we compute the quartet and doublet PES for the N atom 
approaching the Si02 network. Our results are shown in 
Fig.|31 From Eq.|21we obtain the probability for the spin 
conversion using for V the quartet-to-doublet spin-orbit 
matrix element of the N atom {'^P — > ^P), obtained from 
spectroscopy data of 13.8 cm^^. Thus, we find a quartct- 
to-doublet conversion probability Pqd = 1-6 x 10^^. This 
probability corresponds to a single passage of the N atom 
through the crossing point. To estimate the number of 
spin-conversion events of thermally activated N atoms 
approaching the Si02 network, we calculate the ratio be- 
tween the rate of N following the quartet diffusion PES 
(Tq) and the rate of N changing from doublet to sin- 
glet PES at the crossing point (Tqd). To calculate F, 
we use Eqs. 13 substituting AEt by the energy barrier 
for the quartet diffusion PES obtained from Fig. |2Ib) 
{AEq = 0.6 eV). Tqd was estimated by Eqs. 0] substitut- 
ing Pts by Pqd and using AEp = 0.6 cV, the difference 
between the energies at the interstitial site and at the 
quartet-doublet crossing point, obtained from Fig. [3 For 
a temperature of 800 °C, we find that Tq/Tqd ~ 6 x 10^. 
This results indicates that the N atom will most likely 
proceed along the quartet PES, i.e., it diffuses through 
the Si02 essentially without reacting with the network. 

Although very unlikely, the reaction of a doublet N 
with the oxide due to the spin exchange may also oc- 
cur. Again, following the same procedure described 
for NH, we have three possibilities for the doublet N 
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difFusion-reaction mechanism once it is incorporated into 
the oxide: (i) The incorporated N diffuses through the 
Si02 network by a hopping mechanism. From Eq. |3 
substituting AEt by the a calculated hopping energy 
barrier {AEh — 1-8 eV), we find a reaction rate of 
F/i =ro(3.5xl0^^). (ii) N is released from the net- 
work following the doublet PES. From Eq. |S1 substituting 
AEs by the the doublet PES energy barrier, obtained 
from Fig. m {AEd^O.9 eV), and Pt^ by Pqa, we find 
Fj; = Fo(5.9xlO^^). (iii) N is released from the network 
changing from the doublet to the quartet PES at the cross- 
ing point. From Eq. ^ substituting Pts by Pqd and us- 
ing AEp = 0.6 eV, the difference between the energies at 
the interstitial site and at the quartet-doublet crossing 
point, obtained from Fig.El we find F^g = ro(2.4x 10~^). 
Therefore, if a N atom is incorporated into the oxide 
by a quartet-to-doublet spin exchange, a rare event, it 
will be released from the network following the singlet 
PES, according to the reaction (ii). However, as the in- 
corporation reaction of the doublet N is exothermic, it 
will be reincorporated in the network if the timescale 
for relaxation back to the quartet ground state is long 
enough. Finally, we find that the other two possible re- 
actions are very unlikely as can be verified by the ratios 
Td/Th ~ 1.7 X 10^ and Td/Tdq ~ 2.5 x 10^. 

IV. CONCLUSIONS 

We conclude that because of the ionic character of 
Si02, charged species are trapped by the network, where 
positive (negative) species tend to attach to oxygen (sil- 
icon) atoms, forming strong bonds with binding energies 
as high as 3 eV. This suggests a mechanism for the N 
incorporation in bulk Si02 observed in plasma assisted 
method, which would be associated to charged species 
surviving the first stages of the incorporation reactions 
occurring at the surface. Of course, these charged species 
may also transfer the charge to the network becoming 
neutral. According to our results, neutral species do 
not react with the network remaining at interstitial sites, 



therefore, they would be diffusing species in Si02 able to 
reach the Si/Si02 interface. 

We estimate the diffusivities of neutral species through 
a-quartz. Our results show that NH and atomic N have 
increasing diffusivities with temperature, and the high- 
est one among the nitriding species, suggesting that they 
would be efficient for near interface nitridation. This re- 
sult is in good agreement with the finding of high N con- 
centration at the Si/Si02 interface obtained by thermal 
processing in NH3. On the other hand, NO and N2 show 
relatively constant diffusivities in a-quartz for the same 
range of temperatures. This result is also in good agree- 
ment with the finding of relatively low N concentration 
at the Si/Si02 interface obtained by thermal processing 
in N2O, which is the main source of NO and N2 species. 

We also study the incorporation reaction of NH and 
atomic N with the Si02 network driven by a spin ex- 
change mechanism. We find that, for a typical anneal- 
ing temperature (T = 800 °C), a NH molecule diffusing 
through the triplet PES will have, on average, one in 
forteen possibilities to be incorporated into the network 
by a triplet-to-singlet conversion. The incorporated sin- 
glet NH will most likely return to the triplet PES after 
an inverse spin conversion. We also find that atomic N 
may be incorporated into the network by a quartet-to- 
doublet conversion while diffusing through the quartet 
PES, however, only an average of one in 6x10^ events 
will be successful. If this rare event occurs, the incorpo- 
rated N will be released from the network following the 
doublet PES and, subsequently, reincorporated into the 
network due to the doublet exothermic process. How- 
ever, this mechanism will be valid if the timescale for the 
relaxation to the quartet ground state is long enough. 
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